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Abstract 



We present a SUSY SU(5) x T" unified flavour model with type I see-saw mechanism of 
r*j ' neutrino mass generation, which predicts the reactor neutrino angle to be #13 ~ 0.14 close 

p |- to the recent results from the Daya Bay and RENO experiments. The model predicts also 

values of the solar and atmospheric neutrino mixing angles, which are compatible with the 
, existing data. The T 1 breaking leads to tri-bimaximal mixing in the neutrino sector, which 

,JZj ' is perturbed by sizeable corrections from the charged lepton sector. The model exhibits 

geometrical CP violation, where all complex phases have their origin from the complex 
' Clcbsch-Gordan coefficients of T". The values of the Dirac and Majorana CP violating phases 

^ , are predicted. For the Dirac phase in the standard parametrisation of the neutrino mixing 

matrix we get a value close to 90°: 5 = ir/2 — 0.45 9 C = 84.3°, 9 C being the Cabibbo angle. 
The neutrino mass spectrum can be with normal ordering (2 cases) or inverted ordering. In 
each case the values of the three light neutrino masses are predicted with relatively small 
uncertainties, which allows to get also unambiguous predictions for the neutrino- less double 
beta decay effective Majorana mass. 
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1 Introduction 



Understanding the origin of the patterns of neutrino masses and mixing, emerging from the 
neutrino oscillation, 3 H (3— decay, etc. data is one of the most challenging problems in neutrino 
physics. It is part of the more general fundamental problem in particle physics of understanding 
the origins of flavour, i.e., of the patterns of the quark, charged lepton and neutrino masses and 
of the quark and lepton mixing. 

At present we have compelling evidence for existence of mixing of three light massive neutri- 
nos Vi, i = 1,2, 3, in the weak charged lepton current (see, e.g., [T]). The masses of the three 
light neutrinos Ui do not exceed approximately 1 eV, rrii ^ 1 eV, i.e., they are much smaller 
than the masses of the charged leptons and quarks. The three light neutrino mixing is described 
(to a good approximation) by the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) 3x3 unitary 
mixing matrix, C/pmns- In the widely used standard parametrisation [T], C/pmns is expressed in 
terms of the solar, atmospheric and reactor neutrino mixing angles #12, #23 and #13, respectively, 
and one Dirac - 5, and two Majorana [2] - f3\ and /?2, CP violating phases: 

C/PMNS = U= V(9 12 , 023, #13, 5) Q{p x ,h) , (1.1) 

where 

(I \ / C13 s 13 e' iS \ / c 12 S12 \ 

V = c 23 s 23 1 -sn cia , (1.2) 

\ -S23 C 23 / V S l3 e iS C13 / V 1/ 
and we have used the standard notation Cij = cos 9ij , Sij = sin 9ij and 

Q = Diag(e- i/3l/2 ,e- i/32/2 ,l). (1.3) 

The neutrino oscillation data, accumulated over many years, allowed to determine the pa- 
rameters which drive the solar and atmospheric neutrino oscillations, Attt-q = Am^, #12 and 
I Am 2 ^ I = l^ m iil — l^ m 32l, ^23, with a rather high precision (see, e.g., [T]). Furthermore, there 
were spectacular developments in the last year in what concerns the angle 613. In June of 2011 
the T2K collaboration reported [3] evidence at 2.5a for a non-zero value of #13. Subsequently 
the MINOS [1] and Double Chooz [5] collaborations also reported evidence for #13 ^ 0, although 
with a smaller statistical significance. Global analysis of the neutrino oscillation data, includ- 
ing the data from the T2K and MINOS experiments, performed in 0, showed that actually 
sin #13 7^ at > 3a. In March of 2012 the first data of the Daya Bay reactor antineutrino 
experiment on 613 were published [7j. The value of sin 2 2^13 was measured with a rather high 
precision and was found to be different from zero at 5.2a: 

sin 2 20i3 = 0.092 ± 0.016 ± 0.005 , 0.04 < sin 2 20 13 < 0.14 , 3a , (1.4) 

where we have given also the 3a interval of allowed values of sin 2 2#i3 . Subsequently, the RENO 
experiment reported a 4.9a evidence for a non-zero value of #13 [8], compatibe with the Day Bay 
result: 

sin 2 2#i3 = 0.113 ± 0.013 ± 0.019 . (1.5) 

1 This parametrization differs from the standard one. We use it for "technical" reasons related to the fitting 
code we will employ. Obviously, the standard one can be obtained as Diag(l, e 1 0121 , e iQ31 ) = e i/3l/2 Q, with 
«2i = /3i - P2 and Q31 = Pi ■ 
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Parameter 



best-fit (±ler) 



Am| [10- 5 eV 2 ] 
\Am\\ [10" 3 eV 2 ] 



7.62±0.19 7.12 - 8.20 

2.53±g;?g 2.26 - 2.77 

-(2.40tao?) -(2.15-2.68) 

0.320±{j;[j^ 0.27 - 0.37 

°- 49 iao2 0.39-0.64 




sin 



2 




6- 



12 



sin 



2 



'23 



Sill 



2 







13 



0.026lg;gg| 0.015-0.036 
0.027±g:oot 0.016-0.037 



Table 1: The best-fit values and 3<r allowed ranges of the 3-neutrino oscillation parameters 
derived from a global fit of the current neutrino oscillation data, including the Daya Bay and 
RENO results (from [TO]). These values are obtained using the "new" [TT] reactor v e fluxes. If 
two values are given the first one corresponds to normal hierarchy and the second one to inverted 
hierarchy. 

The results on #13 described above will have far reaching implications for the program of future 
research in neutrino physics (see, e.g., [5]). 

A recent global analysis of the current neutrino oscillation data, in which the Daya Bay and 
RENO results on #13 are also included, was published [10] ■ In Table Q] we show the best fit 
values and the 99.73% CL allowed ranges of Am| l5 sin 2 #12, | Ato 2 -^^ | , sin 2 6*23 and sin 2 #13, 
found in [TO]. 

Stimulated by the fact that all three angles in the PMNS matrix are determined with a 
relatively high precision, we report in the present article an attempt to construct a unified 
model of flavour, which describes correctly the quark and charged lepton masses, the mixing 
and CP violation in the quark sector, the mixing in the lepton sector, including the relatively 
large value of the angle #13, and provides predictions for the light neutrino masses compatible 
with the existing relevant data and constraints. The unified model of flavour we are proposing 
is super symmetric and is based on ££7(5) as gauge group and T' as discrete family symmetry. 
It includes three right-handed (RH) neutrino fields Nm, I = e,/j,,r, which possess a Majorana 
mass term. The light neutrino masses are generated by the type I see-saw mechanism [12] and 
are naturally small. The corresponding Majorana mass term of the left-handed flavour neutrino 
fields vil(x), I = e, /i,r, is diagonalised by a unitary matrix which, up to a diagonal phase 
matrix, is of the tri-bimaximal form |13| : 



In order to account for the current data on the neutrino mixing, and more specifically, for the 
fact that #13 7^ 0, £7 TB m has to be "corrected". The requisite correction is provided by the 
unitary matrix originating from the diagonalisation of the charged lepton mass matrix M e (for a 
general discussion of such corrections see, e.g., [14H16] ). Since the model is based on the ££7(5) 
GUT symmetry, the charged lepton mass matrix M e is related to the down-quark mass matrix 
Md- As a consequence, in particular, of the connection between M e and Mj, the smallest angle 
in the neutrino mixing matrix #13, is related to the Cabibbo angle 9 C : sin 2 6q 3 = C 2 (sin 2 9 c )/2 = 




v/V3 




(1.6) 



2 



(sin 2 e c )/2.5, where C ^ 0.9 is a constant determined from the fit. 

The down-quark mass matrix Mj, and the charged lepton mass matrix M e , by construction 
are neither diagonal nor CP conserving. The matrix M e is the only source of CP violation in 
the lepton sector. Actually, the CP violation predicted by the model in the quark and lepton 
sectors is entirely geometrical in origin. This aspect of the SU(5) x T' model we propose is 
a consequence, in particular, of one of the special properties of the group namely, that 

its group theoretical Clebsch-Gordan (CG) coefficients are intrinsically complex }20| . The idea 
to use the complexity of the Clebsch-Gordan (CG) coefficients of T' to generate the requisite 
CP violation in the quark sector and a related CP violation in the lepton sector was pioneered 
in |21j . For the class of models where the CP violation is geometrical in origin, it is essential 
to provide a solution to the vacuum alignment problem for which all the flavon vevs are real. 
In this paper we present a solution of this problem for the models based on the SU(5) x T' 
symmetry. 

Let us note finally that a model of flavour based on the symmetry group SU(5) x T' was 
proposed, to our knowledge, first in [22] and its properties were further ellaborated in |21j 
and [23]. Although some generic features, as like the connection between the reactor mixing angle 
#13 and the Cabibbo angle 9°, which are based on the underlying SU(5) symmetry, are present 
both in the model constructed in |2ip22j and in the model presented here, the detailed structure 
and the quantitative predictions of the two models are very different. The quark, charged lepton, 
RH neutrino mass matrices and the matrix of the neutrino Yukawa couplings have different 
forms in the two models. This leads to considerable differences in the predictions for various 
observables. In the quark sector, for instance, the value of the CKM phase we find is in much 
better agreement with experimental data. More importantly, in the model proposed in [21LI22] . 
the reactor mixing angle #13 is predicted to have the value sin #13 = sin c /(3\/2) = 0.016, which 
is ruled out by the current data on #13. In contrast, due to non-standard SU(5) Clebsch-Gordan 
relations between the down- type quark and the charged lepton Yukawa couplings |16U24j , we get 
a realistic value for this angle. Moreover, in the model we propose both neutrino mass spectra 
with normal and inverted ordering are possible, while the model developed in [21|,l22j admits 
only neutrino mass spectrum with normal ordering [23] . 

The paper is organized as follows. Section 2 is a brief overview of the considered model. In 
section 3 we discuss the quark and charged lepton sector including a \ 2 fit to the experimental 
data. Section 4 is completely devoted to the neutrino sector. There we describe in detail the 
predictions for the mixing parameters (including CP violating phases), the mass spectra and 
observables such as the sum of the neutrino masses, the neutrinoless double beta ((/3/3)o^-) 
decay effective Majorana mass and the rephasing invariant related to the Dirac phase in the 
PMNS matrix, Jcp- We summarize and conclude in section 5. In the Appendix we discuss the 
properties of the discrete group T', the messenger sector which generates the effective operators 
for the Yukawa couplings, and the superpotential, solving the flavon vacuum alignment problem. 

2 Matter, Higgs and Flavon Field Content of the Model 

In this section we describe the matter, the Higgs and the flavon content of our SU(5) x T' 
unified model of flavour. A rather large shaping symmetry, Z12 x Zf x Z| x Z4, is needed to 
solve the vacuum alignment issue and forbids unwanted terms and couplings in the superpotential 
(specifically in the renormalisable one as described in Appendix[Bj as well as in the effective one 

2 There have been also T' models without a GUT embedding, e.g. [17H19] . 
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Table 2: Matter and Higgs field content of the model including quantum numbers. 



after integrating out heavy messenger fields). We further impose an additional U(1)r symmetry, 
the continuous generalisation of the usual imparity. The messenger fields and auxiliary flavons 
used for the flavon superpotential are discussed in the Appendix. 

The model includes the three generations of matter fields in the usual 5 and 10, represen- 
tations of SU(5), F = (d c ,L)L and T = (q,u c ,e c )L and three heavy right-handed Majorana 
neutrino fields N, singlets under SU(5). The light active neutrino masses are generated through 
the type I seesaw mechanism |12j . Furthermore we introduce a number of copies of Higgs fields 
in the 5 and 5 representation of SU(5) which contain as linear combinations the two Higgs 
doublets of the MSSM. To get realistic mass ratios between down-type quarks and charged lep- 
tons [23] and to get a large reactor mixing angle [16] we have introduced Higgs fields in the 
adjoint representation of SU(5) which are as well responsible for breaking the GUT group. 

The matter and Higgs fields including their transformation properties under all imposed 
symmetries are summarised in Tab. [2j Note that the right-handed neutrinos N and the five- 
dimensional matter representations are organised in T' triplets, while the ten-plets are organised 
in a doublet and a singlet. On the one hand this will give us tri-bimaximal mixing (TBM) in the 
neutrino sector before considering corrections from the charged lepton sector and on the other 
hand the complex Clebsch-Gordan coefficients for the doublets will give us CP violation in the 
quark and in the lepton sector finally. 

There are 13 flavons, which will give us the desired structure for the Yukawa couplings that 
will be discussed in the next section. First of all we have three triplets which will develop vevs 
into two different directions in flavour space, 



(0) = U , (0) = Uo , (0 = 1 Uo • (2.1 






The first two flavons will be relevant for the quark and the charged lepton sector and the third 
one couples only to the neutrino sector. 

Then we have introduced four complex T' doublets. Notice that this spinorial representations 
of the T' group are essential since, having complex Clebsch-Gordan coefficients (see AppendixIA")). 
it is responsible of the CP violation in both quark and charged lepton sector. We assume that 
CP is conserved on the fundamental level (all couplings are real) and all flavon vevs are real. 
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Table 3: Flavon fields coupling to the matter sector including their quantum numbers. In fact, 
C' does not couple directly to the matter fields, but it behaves very similar like the other flavons 
and not like the auxiliary flavons e which will be introduced in Appendix [Cj 

In Appendix ICl we give a superpotential that has the desired flavon vev directions as a solution 
and also fixes the phases of the vevs up to a few discrete choices. For the doublets we find the 
vev alignments 



Furthermore we have introduced six flavons in one-dimensional representations of T' which 
receive all non- vanishing (and real) vevs 



(C')=C'o, (0=CSf, <0 = Co, <C">=Cd, ( P )=Po, (P) = P0- (2.3) 



All flavons including their quantum numbers are summarised in Tab. [3l As we will see soon the 
flavon field does not directly couple to the matter sector. Nevertheless, we mention it here 
because it behaves differently than the auxiliary e flavons which we have introduced to get the 
desired alignment and make all vevs real, see Appendix [Ul 

3 The Quark and Charged Lepton Sector 

In this section we describe the superpotential of the quark and charged lepton content of the 
chiral superfields of the model under study. We will consider the three generations of matter 
fields in the usual 5 and 10, five and ten- dimensional, representations of SU(5), F = (d c ,L)i 
and T = (q,u c ,e c )i- The elements of the Yukawa coupling matrices are generated dynamically 
through a number of effective operators which structure is tightly related to the matter fields 
assignment under the T' discrete symmetry. Indeed the Yukawa coupling matrices can be written 
only after the breaking of the T' discrete symmetry. As will be clear soon, in this description CP 
violation in the quark and charged lepton sector is entirely due to geometrical origin, specifically 
from the use of the spinorial representation of the T' group. Finally, in this section we will 




(2.2) 
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present a \ 2 fit analysis that has been performed by us to get the low energy masses and mixing 
parameters in the quark and charged lepton sector. We show as well that the the simple CKM 
phase sum rule from [25] can be applied here. 



3.1 Effective Operators and Yukawa Matrices 

Before we come to the effective operators which will give us the Yukawa couplings we first 
fix the conventions used for the Yukawa matrices. Throughout this paper we will use the RL 
convention, i.e., 

-£ = Y tJ J%f L H + R.c. (3.1) 

or in other words we have to diagonalise the combination Y^Y . Keep also in mind that F = 
(d c ,L) L 8xidT=(q,u c ,e c ) L . 

We restrict ourselves to effective operators up to mass dimension seven. These operators 
generate Yukawa couplings of the order of 10~ 5 or smaller (see our fit results in Tab. [5j). Higher 
dimensional operators hence can be expected to give only negligible corrections. 

After integrating out the heavy messenger fields, see Appendix [0 we obtain the effective 
operators 

+ ^(Tj) 2 ,(H^C')r(^(Tj>) 3 ) 2 , + ^((T^KVtff(C"(^) 2 ")2' , 

which give the up-type quark Yukawa matrices after the flavons developed their vevs. Here A u 
stands for the messenger scale suppressing the non-renormalisable operators in the up-sector and 
in the down-sector we will introduce A^ correspondingly. We have also given the T' contractions 
as indices on the round brackets. Note that in general there are many different contractions 
possible (for T' and to a less degree for SU(5)) which give different results. Nevertheless, 
we have specified in Appendix [B] the fields mediating the non-renormalisable operators which 
transform in a specific way under T' such that we pick up only the contractions which we want. 

Multiplying the T' and SU(5) indices out we obtain for the up-type quark Yukawa matrix 
at the GUT scale (which is roughly equal to the scale of T 1 breaking) 











ud, 







Y u = \ ib u c u ud u , (3.3) 



where ui = (1 + i)/v / 2 and u = (1 — i)/y/2. The parameters a u , b u , c u , d u and e u are (real) 
functions of the underlying parameters. Note at this point, that the phases of the flavon vevs 
have to be fixed. Otherwise the coefficients in the Yukawa matrix are complex parameters and 
we would not be able to make definite predictions anymore. 

For the down-type quarks and charged leptons (remember that those two sectors are closely 
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related in SU(5)) we find for the superpotential 



(d) _ (d) 

d d 

(d) _ (d) 

+ ^(((T a ^ 4 ) 2 ,(^ / ) 2 ")3^)2"(^ 3 V')2' + ^((^ / )2"(C"^ 1) )l»C")2(T a 0) 2 (3.4) 
d d 

(d) 

d 

where we have again specified the T' contractions. From this superpotential and considering the 
correct SU(5) contractions, which we could not display here for the sake of readability, we get 
the down-type quark and charged lepton Yukawa matrices 

/ua d ib' d 0\ (-\uia d Ob d \ 

Y d =\u)b d c d and Y e = \6ib' d 6c d , (3.5) 
V dj \ -\d d ) 

where a d , b d , b' d , c d and d d are (real) functions of the underlying parameters. 

Note that the prediction from the minimal SU(5) model Y d = Yj is broken. Indeed it has 
to be broken to get realistic fermion masses. For the second generation this is known for a 
long time [26] , In some recent work [24] some new relations to fix this issue were proposed. 
From those we will use here y T /yb = —3/2 and y^/y s ~ 6 where y T , y^, yb and y s stand for 
the eigenvalues of the Yukawa matrices associated to the masses of the r, the fi, the b and the 
s quark respectively. Furthermore it was shown in [16j (see also [27]) that those new SU(5) 
Clebsch-Gordan coefficients might also give a large reactor neutrino mixing angle #13. For the 
current paper we have chosen one of the possible combinations given in [16] but we remark that 
in principle also other combinations are still possible which might be realised in another unified 
flavour model with a similar good fit to the fermion masses and mixing angles. 



3.2 Fit Results and the CKM Phase Sum Rule 

In the last section we have discussed the structure of the Yukawa matrices in the quark and 
the charged lepton sector. These matrices have five free parameters, which in principle can 
be fitted to the low-energy mass and mixing parameters using the renormalisation group. But 
doing so one has to take into account SUSY threshold corrections [28] which modify the masses 
and mixing angles significantly. For example without including them, the GUT scale Yukawa 
coupling ratio, y T /yb, would be roughly 1.3 which is not close to the usual GUT prediction 
of 1. There is a large amount of literature on how to use SUSY threshold corrections to get 
b — r Yukawa unification, for recent papers see, for instance, |24 ]l29ll3"0] . From these studies 
it is known that in order to get b — r Yukawa unification, it is necessary to either consider a 
negative //-term or to have a very high - O(10 TeV), SUSY scale. Nevertheless, we will not use 
unifcation but instead we use the recently proposed GUT scale relation y T /yb = 3/2 induced by 
the vev of an adjoint of SU (5) [24] . which is viable in a large region of the parameter space even 
in constrained MSSM scenarios. 

Due to the importance of the threshold corrections for our fit we briefly revise the most 
important formulas which also defines our parameterisation. In [31] the approximate matching 
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Parameter 


Value 


a„ 


5.81 • 10" 6 




-9.96 • 10~ 5 


('a 


-8.55 • 10~ 4 


d„ 


1.99 • 1(T 2 


6 u 


0.525 


(Id 


-2.82 • 10" 5 


b d 


-5.73 • 10" 4 




-5.09 • 10" 4 




2.50 • 10" 3 


dd 


1.82 • 10" 1 


e q tan /3 


0.1788 


tA tan j3 


-0.0001 



Table 4: Values of the effective parameters of the quark and charged lepton Yukawa matrices 
for tan/3 = 35 and -Msusy = 750 GcV. The two parameters e q and parameterise the SUSY 
threshold corrections. The numerical values are determined from a x 2 -fit to experimental data 
with a lowest x 2 per degree of freedom of 2.76. 



conditions at the SUSY scale, Msusy , 



^ T = (l + e /tan^)^ ss T M cos/3, (3.6) 
4^Hl + ^tan/3) yd M f M cos/3, (3.7) 
y! M = (1 + (e, + e A ) tan/3) y 6 MSSM cos/3 , (3.8) 



for the Yukawa couplings and 



l + e 9 tan/3 ^mssm 



13 l + (eg + ej4 )tan/3 l3 y ' 

nSM _ zjMSSM 
7 12 — "12 



(3.10) 

rSM _ A MSSM fo i i \ 

°CKM — °CKM i K 6 - 11 ) 

for the quark mixing parameters were given, where the SUSY threshold corrections are 
parameterised in terms of the three parameters ei, e q and e^. We will adopt this parameterisation 
neglecting ei, which is usually one order smaller than e q [29]. Furthermore we want to assume, 
that SUSY is broken similar to the constrained MSSM scenario with a positive /i parameter and 
hence we adopt the recently proposed GUT relation y T /yb = 3/2 for the third generation, as 
mentioned earlier. For the second generation we use y^/ys ~ 6 |24j . 

We have fixed the SUSY scale to 750 GeV, the GUT scale to 2 x 10 16 GeV and tan p to 35. 
Therefore we have to fit the ten parameters in the Yukawa matrices and the two parameters 
from the SUSY threshold corrections to the thirteen low energy observables in the quark and 
the charged lepton sector (nine masses, three mixing angles and one phase), so that we have one 
prediction (degree of freedom). 

The RGE running and diagonalisation of the matrices was done using the REAP package [32] . 
Performing a x 2 fit we have found as minimum the results listed in Tab. H] for the parameters 
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Quantity (at mt(mt)) Experiment Model Deviation 



y T in 1(T 2 1.00 0.99 -0.388 

y M in 10" 4 5.89 5.90 0.044 

y e in 10- 6 2.79 2.79 -0.003 



Vb in 10' 
Vs in 10' 

Vs/Vd 


-2 
-1 


1.58 ±0.05 
2.99 ± 0.86 
18.9 ±0.8 


1.57 
2.57 
18.9 


-0.157 
-0.484 
-0.012 


yt 

y c in 10" 
y u in 10 


-3 
-6 


0.936 ± 0.016 
3.39 ± 0.46 
7 01+ 2 - 76 


0.936 
2.79 
7.01 


0.0001 
-1.317 
-0.0003, 


0CKM 

/iCKM 
^13 

ScKM 




n OOK7+0.0009 
_ .0010 

041 ^+ - 0011 
0.0036 ± 0.0002 

1 on9q+0-0786 
±.^U^O_g 0431 


0.2257 
0.0416 
0.0036 
1.2610 


-0.0107 
0.1268 
0.2043 
0.7465 



Table 5: Fit results for the quark Yukawa couplings and mixing and the charged lepton Yukawa 
couplings at low energy compared to experimental data. The values for the Yukawa couplings are 
extracted from [33], the ratio y s /yd is taken from [34] and the CKM parameters from pQ. Note 
that the experimental uncertainty on the charged lepton Yukawa couplings are negligible small 
and we have assumed a relative uncertainty of 3 % for them. The x 2 per degree of freedom is 
2.76. A pictorial representation of the agreement between our fit and experiment can be found 
as well in Fig. [T] 



X 2 /dof = 2.76 



+2(T 

+ \<r 



Oo- • • 



* " . 



• • 



. . T 



-lcr 
-2a 



e fi T — S t> U C 1 * 12 "l3 "23 "CKM 



Figure 1: Pictorial representation of the deviation of our fit from low energy experimental data 
for the charged lepton Yukawa couplings and quark Yukawa couplings and mixing parameters. 
The deviations of the charged lepton masses are given in 3% while all other deviations are given 
in units of standard deviations a. 
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and in Tab. [5] and in Fig. [I] we have presented the results of the fit for the low energy observables 
compared to the experimental results. Note that we have assumed an uncertainty of 3% on the 
Yukawa couplings for the charged leptons. Their experimental uncertainty is much smaller, so 
that their theoretical uncertainty (Accuracy of RGEs, neglecting SUSY threshold corrections 
for the leptons, NLO effects, ...) is much bigger, which we estimate to be 3%. 

We find good agreement between our model and experimental data with a minimal x 2 P er 
degree of freedom of 2.76. In fact this agreement is not accidental. We have chosen the SU (5) 
coefficients such that, we expect good agreement and we have also enough free parameters to fix 
the mixing angles. In other words one could determine the eigenvalues and mixing angles from 
the data and then the CKM phase would be a prediction. But as we will demonstrate now, the 
choice for our phases in the Yukawa matrices was done in such a way, that we can expect a good 
prediction for the CKM phase as well. 

We will show in the following that the sum rule given in [25] can be used here. To apply the 
sum rule we have to find approximate expressions for the complex mixing angles (sec [25]). For 
the rest of the subsection we will use the the notation of [25] which we just briefly summarise 
here for convenience. The CKM matrix Uqkm can be written as 



f^CKM = U Ur U 



L^d L ~ \ u 23 u 13 U 12 I u 23 u 13 u 12 ' 



(3.12) 



where the matrices U UL and Ud L diagonalise the up- and down-type quark mass matrices and 
the unitary matrix 

-1612 qN 



U 



12 



COS 012 

sin #12 e 




i<5i2 



sin 0i2 e 

COS 012 









(3.13) 



The matrices U\% and C/23 are given by analogous expressions. 

We find at leading order for the respective mixing angles and phases 



?f 2 e~ i<5 « 



'12 



C ld 12 



'13 



nd 
; 23 







4 



'23 



da 



'13 



e 



bud-u 



(3.14) 
(3.15) 



where we have used for 0" 2 that d\ 



So we see that 5f 2 is 



l/2c u and e u ~ 0.5 from our fit. 
not simply tt/2 as one would expect from a quick first inspection. Note also that the phase sum 
rule was derived for 0" 3 = 6f 3 = 0, which is not exactly true in our case for 0" 3 . But in fact it 
is sufficient, that <C #i 2 #23 which is fulfilled here. 

The angle a in the CKM unitarity triangle is experimentally measured to be a = (90.7 
1351 for which the sum rule 

^2 , (3-16) 



7+4.5 > 
-2.9J 



a 



°12 



was given in |25J . Plugging in our approximate analytical expressions for , eqs. (|3. 14)) and 
(|3.15p . we find that a w tt/2 and our model is in good agreement with experimental data as we 
have also seen it before from our numerical fit. 



4 Neutrino Sector 



The model includes three heavy right-handed Majorana neutrino fields iV which are singlets 
under SU (5) and a triplet under T'. Through the type I seesaw mechanism [12] we generate light 
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neutrino masses. The neutrino sector is described by the following terms in the superpotential 



W l/ = X 1 NN^ + NN(X 2 p + X 3 p) + ^-(NF) 1 (HPp) 1 + ^-(NF) 1 (HPp) 1 , (4.1) 

where we have given the T' contractions as indices at the brackets for non-renormalisable terms 
and from now on A labels a generic messenger scale. Note that the contraction of three triplets in 
general is not unique, see also Tab. because the product of two triplets contains a symmetric 
and an antisymmetric triplet. But since we multiply here two N with each other only the 
symmetric combination gives a non-vanishing contribution. In the following we will discuss the 
phenomenological implications of this superpotential (including corrections from the charged 
lepton sector). 



4.1 The Neutrino Mass Spectrum 

From eq. (|4.1|) we obtain for the mass matrix for the right-handed neutrinos and the Dirac 
neutrino mass matrix 

-Z -Z \ A 0\ ( 

2Z -Z + l] , M c =|0 




-z + x 



where X, Z and p' are real parameters depending on the couplings and the vevs in eq. (|4.ip . 
The right-handed neutrino mass matrix Mr is diagonalised by the tri-bimaximal mixing (TBM) 
matrix [13] 

/ ^/2/3 y/T/3 \ 
U TBM = \-y/l/6 y^T/3 -7T72 , (4.3) 

V-x/iTe 7173 v/172/ 

such that the heavy RH neutrino masses read: 

U^ BM M R U TBM = D N = Di ag (3Z+X,X,3Z-X) = Diag(M 1 e i ^ 1 , M 2 e^ 2 , M 3 e i<fe ) , M lj2>3 > 0, 

(4.4) 

where 

Mi = \X + 3Z\ = \X\\l + aj% 0i = arg(X + 3Z) (4.5) 
M 2 = \X\, 2 = arg(X) (4.6) 
M 3 = \X-3Z\ = \X\\l-ctj% 3 = arg(3Z - X) . (4.7) 

Here a = \3Z/X\ > and = arg(Z) — arg(X). Since X and Z are real parameters, the phases 
01) 02) 03 and take values or ir. A light neutrino Majorana mass term is generated after 
electroweak symmetry breaking via the type I see-saw mechanism: 

M v = -Ml M^ 1 M D = U* u Diag (m 1 ,m 2 ,m 3 ) U$ , (4.8) 

where 

f/^i^BMDiag^^/ 2 ,^^ 2 ,^^ 2 ) =iU TBM Q, Q = Diag(e i ^/ 2 ,e i ^ 2 ,e i ^/ 2 ) , (4.9) 
and w.1^,3 > are the light neutrino masses, 

m=K) 2 ±, -1,2,3. (4.io) 
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The phase factor i in eq. (j4.9j) corresponds to an unphysical phase and we will drop it in what 
follows. Note also that one of the phases fa, say fa, is physically irrelevant since it can be 
considered as a common phase of the neutrino mixing matrix. In the following we always set 
fa = 0. This corresponds to the choice (X + 3Z) > 0. 

The type of the neutrino mass spectrum in the model is determined d by the value of the 
phase (ft. Indeed, as it is not difficult to show, we have: 



X\ 2 \AJ ll + ae^ril 



Thus, for cos (ft = +1, we get Am^ > 0, i.e., a neutrino mass spectrum with normal ordering 
(NO), while for cos (ft = —1 one has Am^ < 0, i.e., neutrino mass spectrum with inverted 
ordering (10). We have also: 



Am 2 21 = Aml = T ^[! T ) ^ ' ■ ( 4 - 12 ) 



1 / p'\ 4 a (a + 2 cos < 
= Txf { l) |i + ae i 



For a given type of neutrino mass spectrum, i.e., for a fixed (ft = or ir, a constraint on the 
parameter a can be obtained from the requirement that Am^ > and from the data on the 
ratio: 

Am? 1 

r = — — f- = - (a + 2 cos <ft) (l - 2a cos (ft + a 2 ) = 0.032 ± 0.006 . (4.13) 
|Amjy 4 

Using the values of a thus found and the value of, e.g., Am^, one can get (for a given type of 
the spectrum) the value of the factor in eq. (|4.12|) . \X\~ 2 (p' /A) 4 . Knowing this factor and a, 
one can obtain the value of the lightest neutrino mass, which together with the data on Am^ 
and A?ng 1 ^ 32 ^ allows to obtain the values of the other two light neutrino masses. Knowing the 
latter one can find also the two ratios of the heavy Majorana neutrino masses. 

In the case of NO neutrino mass spectrum (0 = 0), there are two values of a which satisfy 
equation (j4TT3j) for r = 0.032: a = 1.20 (solution A), and a = 0.79 (solution B). In the case of 
solution A, as it is not difficult to show, the phases 

fa = , fa = , solution A (NO) , (4.14) 

and the three neutrino masses have the values: 

mi ^ 4.44 x 10" 3 eV ,m 2 = 9.77 x 10~ 3 eV , m 3 ^ 4.89 x 10~ 2 eV , solution A (NO) . (4.15) 

Evidently, the spectrum is mildly hierarchical. The ratios of the heavy Majorana neutrino masses 
read: Mi/M 3 ^ 11.0 and M 2 /M 3 ^ 5.0. Thus, we have M 3 < M 2 < M 1 . 
For solution B we find 

fa = , fa = 7T , solution B (NO) , (4.16) 
while for the values of the three neutrino masses we get: 
mi =■ 5.89 x 10~ 3 eV , m 2 = 1.05 x 10~ 2 eV , m 3 ^ 4.90 x 10" 2 eV , solution B (NO) . (4.17) 



3 We are following in this part the similar analysis performed in |36] 
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Figure 2: The values of the three light neutrino masses corresponding to the solutions A (left 
panel) and B (right panel) in the case of NO spectrum, versus r. The dotted, dashed and solid 
lines correspond to the three light neutrino masses mi, m2, 7713. The gray region is excluded by 
present oscillation data. The vertical dashed line corresponds to the best fit value for r ~ 0.032. 
See text for further details. 



The heavy Majorana neutrino mass ratios are given by: M1/M3 = 8.33 and M2/M3 = 4.67. 
Therefore also in this case we have M3 < M2 < M\. 

For the IO spectrum (<j) = 7r), we find only one value of a which satisfies eq. (|4.13p with 
r = 0.032: a = 2.014. The phases 4>2 and 4>3 take the values: <p2 = tt, 03 = 0. The light neutrino 
masses read: 

mi 9* 5.17 x 10~ 2 eV,m 2 = 5.24 x 10~ 2 eV,m 3 9* 1.74 x 10~ 2 eV , (IO) , (4.18) 

i.e., the light neutrino mass spectrum is not hierarchical exhibiting only partial hierarchy. For 
the heavy Majorana neutrino mass ratios we obtain: M1/M2 = 1.014 and M3/M2 = 3.01. Thus, 
in this case N\ and N2 are quasi-degenerate in mass: M\ = M2 < M3. 

In the Figs. [2] and Owe present the dependence of the neutrino masses with respect to r for 
normal and inverted ordering respectively. 

4.2 The Mixing Angles and the Dirac and Majorana CP Violation Phase 

The PMNS neutrino mixing matrix received contributions from the diagonalisation of the neu- 
trino Majorana mass matrix M u and of the charged lepton mass matrix M e = VdY e : C/pmns = 
ul L U Ul where U v is given in eq. (|4.9p with Q = Diag(l, e 1 ^ 2 / 2 , e 1 ^ 3 / 2 ) and the values of the phases 
4>2 and 03 in the cases of NO and IO spectra were specified in the preceding subsection. The 
matrix of charged lepton Yukawa couplings Y e , eq. (j3.5j) . and thus M e , has a block-diagonal form. 
The unitary matrix U e L diagonalises the the Hermitian matrix M\~M e : M.\M e = U e L{M^) 2 U^ eL , 
where Mg = diag(m e , m^, m T ), mi being the mass of the charged lepton /. As a consequence 
of the block-diagonal form of M e , the matrix U e L can be parametrised in terms of one mixing 
angle (0f 2 ) and one phase (ip): U e L = <I> i?i2(#f 2 ), where $ = diag(l,e 1</3 , 1) and 



-^12 (#12) - 




(4.19) 
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Figure 3: The values of the three light neutrino masses in the case of the solution corresponding 
to 10 spectrum, versus r. The dotted, dashed and solid lines correspond to the three light neutrino 
masses to 3 , mi , to 2 . The gray region is excluded by present oscillation data. The vertical dashed 
line corresponds to the best fit value for r = 0.032. 



Due to the SU(5) symmetry of the model, Y d and Y e (and therefore the corresponding down 
quark and charged lepton mass matrices) are expressed in terms of the same parameters. As a 
consequence, the angle Q\ 2 in the model considered is related to the Cabibbo angle 6 C = 0.226. 
Using, for example, the approximate formulas from [16], we find that 



'12 



6 C ^ 0.9 6 C , 



where we have used the values of b' d and b d from Table HI 

Comparing next the expressions on the two sides sides of the equation M<fM e 
we get, in particular: 



(4.20) 



U eL (M^ul 



2 ' (m — m e ) cos 0\ 2 sin 9 



v d 



\b d a d - 36c d b' d 



(4.21) 



Using the fit results in Table S] one can check that the right hand side of the last equation is 
real and positive. Comparing the phases of the two expressions one concludes that 



-7T 



(4.22) 



In the approximation we are using the PMNS matrix is given by: 

C^PMNS = 



2/3cf 2 + y/lTUs^e- 1 ^ 
1/6 



1/3^-^/1/3^-^ 
l/3 S f 2 + y573c e 12 e-^ 
1/3 



Q , (4.23) 



where cf 2 



cos6f 2 , s| 2 = sin#f 2 and Q is the diagonal phase matrix denned in eq. (|4.9p . It 
follows from the above expression for the PMNS matrix that the angle #13 is given approximately 
by 



sin 0i3 



sin 



2qc 



2.5 



0.02 , 0.£ 



(4.24) 



where we took into account the relation in eq. (|4.20p and the value of C = \b' d /bd\ 
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As was shown in, e.g., [16], the phase ip and the Dirac phase 5 in eqs. (|1.1[) - (jl.2p are related 
(at leading order) as follows: 

<5 = 92 + 71-. (4.25) 
Thus, for the Dirac phase we get from (|4.22|) : 



7T 

2 



(4.26) 



Numerically, for 99 = 3vr/2 and s\ 2 = 0.203 (see eq. K20\i ). the PMNS matrix, eq. (|4T23l . 

16 eld S " 

/ 0.804^ 5 - 81 ° 0.577 e - in - 50 ° 0.144 e - i27aooo ° \ 

C/pmns= 0.433e- i67 - 85 ° 0.577^ 78 ' 50 ° - 0.692e~ i27a000 ° Q. (4.27) 
\ -0.408 0.577 0.707 / 

Thus, comparing the absolute values of the elements U e ±, U e 2, and U T z of the PMNS 
matrix in the standard parametrisation, eq. (jl.l|) . and in eq. (|4.27|) . we have: C12C13 = 0.804, 
S12C13 = 0.577, S23C13 = 0.692 and C23C13 = 0.707. Using the predicted value of #13, eq. (j4.24|) . 
these relations allow us to obtain the values of 612 and #23- We note that the tri-bimaximal 
mixing value of the solar neutrino mixing angle #12, which corresponds to sin 2 #12 = 1/3, is 
corrected by a quantity which, as it follows from the general form of such corrections |14H16| . is 
determined by the angle #13 and the Dirac phase 5: 

1 2 [2 

sin 2 #12 = - + ^— sin 6 13 cos 5 (4.28) 

o o 

where 5 is the Dirac phase in the standard parametrisation of the PMNS matrix. As we have 
seen, to leading order 5 = ir/2. The Majorana phases fix, 02 (or «2i and 0:31) are determined, as 
it follows from eqs. (jl.ip and (|4.23p (or (|4.2T|) ) . by the diagonal matrix Q and take CP conserving 
values. Note, however, that the parametrisation of the PMNS matrix in eq. (|4.27|) differs from 
the standard one: it corresponds to one of the several possible parametrisations of the PMNS 
matrix [TS]. Thus, in order to get the values of the Dirac and Majorana phases 5 and (3i, /?2 (or 
«2i 1 ^31 ) j of the standard parametrisation of the PMNS matrix, one has to bring the expressions 
(|4.27p in a form which corresponds to the "standard" one in eq. (jl.ljl . This can be done by 
using the freedom of multiplying the rows of the PMNS matrix with arbitrary phases and by 
shifting some of the common phases of the columns to a diagonal phase matrix P. The results 
for the numerical matrix in eq. (|4,27p is: 

0.804 0.577 0.144e" i84 - 25 ° 

C/pmns = I - 0.433^ ia59 ° 0.577e- i5 ' 75 ° 0.692 \PQ, (4.29) 

0.408e" ilL56 ° -0.577e i5 ' 75 ° 0.707 

where Q = Diag(l, e i<fe / 2 , e i<fe / 2 ) = e i(fe / 2 Diag(e~ i<fe / 2 , e~ i (*s-«a)/2 j 1) an d the new phase ma- 
trix P = Diag(e llL50 ,e _l5 ' 81 ,—1). Now comparing eq. (|4.29|) with eq. (jl.ip we can obtain 
the values of the Dirac and the two Majorana phases of the standard parametrisation of the 
PMNS matrix, predicted by the model. For the Dirac phase we find 5 = 84.3°. Note that the 
Majorana phases and ^2/2 (or a2i/2 and 0:31/2) in the standard parametrisation are not 
CP conserving [23| : due to the matrix P they get CP violating corrections to the CP conserving 
values and it/ '2 or 3tt/2. 

As we have seen, the value of the Dirac phase 5 predicted by the model is close to n/2. This 
implies that the magnitude of the CP violation effects in neutrino oscillations, is also predicted 
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Quantity 


Experiment {2a ranges) 


Model 


sin 2 9 12 


0.275 - 0.342 


0.340 


sin 2 #23 


0.36 - 0.60 


0.490 


sin 2 6»i3 


0.015 - 0.032 


0.020 


5 




84.3° 



Table 6: Numerical results for the neutrino sector. The experimental results are taken from [5] 
apart from the value for #13 which is the DayaBay result [7] . 



to be be relatively large. Indeed, the rephasing invariant associated with the Dirac phase [37] . 
Jcp = Im(L r * 1 C/ /i iC/ e 3L/* 3 ), which determines the magnitude of CP violation effects in neutrino 
oscillations [38], has the following value: 

J C p = 0.0324 . (4.30) 

The values we have obtained for both sin #13 and 5 are in very good agreement with the 
numerical results in Table [6] derived using the REAP package [32] . 

It is possible to derive simple analytic expressions which explain the numerical results ob- 
tained above and quoted in Table [6] Indeed, up to corrections of order (#f 2 ) 2 we have: 



?12 


1 V^/^p N 9 

= arcsm-^= + — (0 12 ) , 


(4.31) 


h 3 


1 9* 

~ T2 12 ' 


(4.32) 






(4.33) 


6 


- 2 ~ 2 12 ' 


(4.34) 


Pi 


= 2vr - 20f 2 + 3 , 


(4.35) 


ft 


= 2vr + 9\ 2 + 03 - 02 , 


(4.36) 



where 9\ 2 — 0.888# c . Note that the expression for 5 is correct up to 0{9\ 2 ) only because it 
appears always with #13 which is of order 9\ 2 itself. Numerically, these approximations give for 

#12 



' e -0.2: 



(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 

As we see, the results obtained using the approximate analytic expressions are in very good 
agreement with those derived in the numerical analysis. 



sin 2 #i2 


= 0.340 , 


sin 2 6>i3 


= 0.020 , 


sin 2 #23 


= 0.490 , 


5 


= 84.3° , 


Pi 


= 337.1° + 03 , 


Ih 


= 11.5° + 03 "02 
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Note that all these relations were derived neglecting RGE corrections. Indeed they are under 
control. For the inverted ordering the RGE corrections can be expected to be largest, because 
there mi and mi are almost equal [32]. We have found numerically with the REAP package |32j 
that the biggest deviation is in 5 which goes down to 81.2°. The Majorana phases run less than 
one degree and also the mixing angles stay well within their two sigma ranges. 

4.3 Predictions for Other Observables in the Neutrino Sector 

We derive in this section the predictions for the sum of the neutrino masses and the effective 
Majorana mass |(m)| in neutrinoless double beta decay (see, e.g., [3D] ) using the standard pa- 
rameterisation of the PMNS mixing matrix as in (|1.2|) and the results on the neutrino masses, 
mixing angles and CPV phases obtained in preceding subsections of this Section. 

In the case of solution A for the NO neutrino mass spectrum we get for the sum of the 
neutrino masses: 

3 

rn k = 6.31 x 1(T 2 eV , solution A (NO) . (4.43) 

k=l 

In this case we have = 03 = (see subsection 4.1) and for the effective Majorana mass we 
obtain using eqs. (|4.15p and (|4.29p : 

3 

\(m)\ = | ^2(U PM m) 2 ek m k\ = 4.90 x 10~ 3 eV , solution A (NO) . (4.44) 
k=l 

The same quantities for solution B of the NO spectrum have the values: 

3 

^2 m k = 6.54 x 1(T 2 eV , solution B (NO) , (4.45) 
fc=i 

and 

|(m)| = 7.95 x 1(T 3 eV , solution B (NO) , (4.46) 

where we have used the fact that for solution B we have 02 = and 03 = ir. As a consequence, 
in particular, of the values of 02,3, the three terms in the expression for \{m)\ essentially add. 
Finally, in the case of 10 spectrum we obtain: 

3 

Y m k = 12.1 x 1(T 2 eV , (10) , (4.47) 

k=l 

and 

|(m)| = 2.17 x 1(T 2 eV , (10) , (4.48) 

We recall that for the 10 spectrum we have 02 — tt and 03 = and there is a partial compensation 
in \ {m)\ between the dominant contributions due to the terms oc m\ and oc m<i- 

5 Summary and Conclusions 

We have presented here the first SU(5) x T' unified model of flavour, which predicts the re- 
actor neutrino mixing angle #13 to be in the range determined by DayaBay [7] and RENO [H] 
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experiments, and all other mixing angles are predicted to have values within the experimental 
uncertainties. It implements a type I seesaw mechanism and from the breaking of the discrete 
family symmetry T" we obtain tri-bimaximal mixing in the neutrino sector. The relatively large 
value of #13 is then generated entirely by corrections coming from the charged lepton sector. This 
is a generic effect in GUTs where Yukawa couplings are related to each other. Here we have used 
recently proposed SU(5) GUT relations |24| between the down-type quark Yukawa matrix and 
the charged lepton Yukawa matrix to get the relatively large prediction for the reactor mixing 
angle #13 along the lines proposed in [161127], 

The corrections to the solar and the atmospheric neutrino mixing angle are under control 
due to the structure of the charged lepton Yukawa matrix and the pattern of the complex CP 
violation phases. The model exhibits a special kind of CP ciolation, the so-called "geometrical" 
CP violation. All parameters and vevs are real and all non-trivial phases are coming from the 
complex Clebsch-Gordan coefficients of T' and are integer multiples of 7r/4. We have given the 
renormalisable superpotential which generates effectively the Yukawa matrices after integrating 
out heavy messenger fields and plugging in the family symmetry breaking flavon vevs which was 
missing so far in the literature for SU (5) x T' models. The flavon vevs point in special directions 
in flavour space and are all real. These results come out as solutions to the flavon alignment 
superpotential we have presented in the Appendix. 

We have shown, in particular, that the phase pattern in the Yukawa matrices actually gives a 
good fit of the quark and charged lepton masses and the CKM parameters at low energies. This 
fit fixes the charged lepton Yukawa matrix completely and since we find tri-bimaximal mixing 
in the neutrino sector itself, we can make predictions for the neutrino masses and all PMNS 
parameters. The angle #13 is predicted to have a value corresponding to sin 2 #13 = 0.8 sin 2 6 c /2 = 
0.02. For the Dirac phase 5 we obtain in the standard parametrisation of the PMNS matrix 
5 = 84.3°. Our model also predicts sin 2 #12 = 0.340 and sin 2 #23 = 0.490. There are three 
different possible solutions for the neutrino masses, two with normal ordering (NO, solutions 
A and B) and one with inverted ordering (IO). All three cases can be tested in experiments 
determining the absolute neutrino mass scale (or the sum of the three neutrino masses), in 
experiments which can measure the solar and atmospheric neutrino mixing angles with a high 
precision, in experiments searching for CP violation in neutrino oscillations and in neutrinoless 
double beta decay experiments. For the sum of three neutrino masses we get (with relatively 
small uncertainties, see Figs. (2) and (3)): J2k=i m k = 6.31 x 10" 2 eV (NO, A); 6.54 x IO" 2 eV 
(NO, B) and 12.1 x 10~ 2 eV (IO). The (/3/3)o^-decay effective Majorana mass for the three 
solutions is also unambiguously predicted: \(m)\ = 4.90 x IO -3 eV (NO, A); 7.95 x 10 -3 eV (NO, 
B); 2.17 x 10 -2 eV (IO). The three solutions differ only in the values of the three neutrino masses 
and of the Majorana phases, so that we make one single prediction for the rephasing invariant 
which determines the magnitude of CP violation effects in neutrino oscillations: Jcp = 0.0324. 
This value of Jcp is relatively large and can be tested in the experiments on CP violation in 
neutrino oscillations. 

In conclusion, with the recent measurement of the last unknown neutrino mixing angle, 
neutrino physics has entered a new era. All angles are determined with a rather good precision, 
constraining flavour models severely. Since #13 turned out to be relatively large, the observation 
of CP violation in the lepton sector might be feasible with data from the running and upcoming 
neutrino oscillation experiments. Explaining the data on leptonic CP violation would pose 
another challenge for flavour models. The model we proposed here is from this point of view 
rather comprehensive combining many ideas which have been proposed elsewhere but have been 
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combined here consistently for the first time. Due to the GUT structure we can fit the quark 
masses and mixing parameters and the charged lepton masses, and using the latter we make 
definite predictions for the neutrino mass spectrum, the leptonic mixing angles and the leptonic 
CP violating phases. Our model is therefore testable in a variety of experiments. We are looking 
forward to the outcome of these tests. 
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A T': The Rules of the Game 

T' is the double-covering group of the tetrahedral symmetry T which is isomorphic to A4, the 
group of the even permutations of four objects. T 1 contains three inequivalent one-dimensional 
representations, called 1, 1' and 1", one three-dimensional, 3 and three two-dimensional repre- 
sentations, 2, 2' and 2" . Two of these representations are real, 1 and 3, one is pseudo-real 2 
and the other four ones are complex. We list in Tab. [7] the relevant tensor products of T". For 
more details on T', see, e.g. [T7] and references therein. 

B Messenger Sector 

In our model we consider non-renormalisable operators. In general the contraction of the SU (5) 
and T' indices may not be unique which is nevertheless essential for our model. Our predictions 
are based on the fact, that only a certain contraction is allowed as we have, for example, indicated 
in eq. (|3.2p for the T" indices. For the connection between the so-called UV completion and 
predictivity of a model see also |41j . Hence we have to specify the so-called messenger fields 
which generate only the desired contractions in the operators after being integrated out in a 
specific order. 

The full list of messenger fields of our model is given in Tab. [HJ Every messenger pair in 
every line receives a mass term in the superpotential, like, for example, M^jS^S^. For the sake 
of brevity we do not write down all mass terms, but it is important to note, that there are no 
mass terms between messengers in different lines allowed. We assume all the messenger masses 
to be above the scale of T 1 and SU (5) breaking, which are closely related in our model as we 
will see in the next section. Many messengers carry SU(5) quantum numbers so that above the 
messenger scale, which we denote by A the gauge coupling becomes quickly non-perturbative, 
so that we are not predictive above this scale. 

After this general remarks we now turn to the superpotential which describes the couplings 
of the various fields to the messengers. We start with the messengers coupling to the matter, 
Higgs and flavon fields. The supergraphs showing these couplings are given in Figs. HKl From 
these diagrams one can read off all the relevant contractions and couplings. Nevertheless, we 
give now the renormalisable superpotential containing the messenger fields. 
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Table 8: Messenger fields used in our model. After integrating out these fields we end up with 
the desired effective operators. For the sake of brevity we do not list all mass terms in the text. 
The messenger pair in every line has a mass term and there are no cross terms allowed. 
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p,p 



Figure 6: The supcrgraphs before integrating out the messengers for the neutrino sector. 

Apart from the messenger mass terms (which we do not write down explicitly) there are no 
terms with one or two fields involving matter, Higgs and flavon fields. For the down-type quark 



diagrams we find (here and in this whole section we do not write down the couplings) 

W T d cn = Fh£ ] TI + 0T§T$„ + H% 4 T 3 r b r (B.l) 

+ T a <t>%, + H' 2 ' 4 r c 2 r c 2 , + ^'r c 2 r d 3 + FH^ri (b.2) 

+ T a H'^t c 2 , + fl£Vns„ + n^ng + n^T c 2 „E^, (b.3) 

+ H>~? 2 , + E a 2ff T e 2 M\ ff + ("H^nl, + ("? e 2 ,? b 2 + T a ^T b 2 (B.4) 

+ Fi;"Z a 2 „ + S^Hj,r§, + H' 2 ^"r a 2 „ + fl£2j„T?„ + T^Taip" , (B.5) 

for the up-type quarks 

W.r 1 = H^ ] Tl + T a ^T a 2 „ + ^ 2) T|,T|„ + r 3 ^"T|, (B.6) 

+ T a 4>"Tl + C'^s 1 ^?" + ni'T«»T<| (B.7) 

+ r <£Tf„ + r§>nj,T|, + flf ^n?* + ^"r| w T| + t^t^' (b.8) 

+ T^T^ + T£C"T^ + ff 5 (3) T^Tf, + C"T^T^, + T a (f>T 2 , , (B.9) 
and for the neutrino sector 

yyrcn = ^ + ^ + ^2 ~ + ^jy-jja + P Q,\ + H^ 2) pClf . (B.10) 

There are five additional operators which generate dimension eight or more operators in the 
matter sector, which we neglect. For completeness we give them as well 

w ren, matter = g-fc-fc, + TgTfV" + V"^2»^3 + f 2'^2"^3 + f . (B.ll) 

We turn now to the messengers, which give the non-renormalisable terms in the flavon 



alignment superpotential which we denote collectively with S. In this sector all the supergraphs 
have the structure as given in Fig. [JJ (The role of the auxiliary e fields is described in the 
next section and their quantum numbers are given in Tab. [TUJ) For the sake of brevity we 



23 




Figure 7: One typical diagram for the messengers in the fiavon sector. We consider only effective 
operators up to dimension four. For the sake of brevity we only show one diagram. The other 
ones are quite similar with the driving held on one side and the auxiliary e helds on the other 
side. 

do not give all the diagrams for the fiavon sector, but all diagrams can easily be derived from 
the renormalisable fiavon superpotential. We give here only the terms, where a messenger is 
involved. The terms, where no messenger is involved will be discussed in the next section, when 
we discuss the superpotential responsible for the fiavon alignment. The superpotential involving 
the £ fields reads 



Apart from these there are as well operators which give dimension five operators in the fiavon 
alignment superpotential after integrating out the messenger fields, which we we will neglect. 
These operators are 



Wd>5 flaV ° n = + S'^l^t, + S£'£f„£?„ + S'lZ\„Tl\„ + S~£g£g + S'fi&l (B.17) 

+ + S'^fri + (S ei + S^ + S p )^ c 3 + (S ei +Ss + S p )^i . (B.18) 



Now we have discussed the messenger sector. After integrating out the messengers from the 
renormalisable superpotential we end up with the effective operators which give us the desired 
fiavon vev alignments and structures of the Yukawa matrices. 

C Fiavon Vacuum Alignment 

In this appendix we present the solution for our fiavon vacuum alignment. In the present model 
all the discussed results crucially depend on the vev structure and on the fact that all fiavon 
vevs are real. 

In the fiavon potential two new kinds of fields are introduced. First we have to add driving 
fields which are gauge singlets but transform in a non trivial way under the family and shaping 
symmetries and have a U(1)r charge of two. Minimising the F-term equations of this fields will 
give us the correct alignment (including phases) as one possible solution. Second we introduce 
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Table 9: List of the driving fields from the superpotential which give the desired vacuum 
alignment. All driving fields are SU(5) gauge singlets and charged under U(1)r with charge +2. 
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Table 10: List of the auxiliary flavon fields that do not couple to the matter sector. The ej 
fields are all SU(5) x T 1 singlets and carry no U(1)r charge. 



auxiliary fields 6j, i = 1, . . . , 13, which are singlets under SU(5) and T', but they transform in a 
non trivial way under the additional shaping symmetries. They appear only in the flavon super- 
potential. Indeed, these fields are introduced to compensate the charges of different operators, 
so that they are related to each other in the F-term equations. Note that we have to include for 
our alignment non-renormalisable operators, where we restrict ourselves to operators with mass 
dimension not higher than four in the superpotential. The driving fields are listed in Table 
and the auxiliary fields are listed in Table [TUJ 

Before going into the more complicated details of the flavon vacuum alignment we briefly 
discuss the "alignment" of the auxiliary flavons, which is simply the question how to give them 
a real vev. For this purpose we used the simple idea advocated in [42] . which we can directly 
illustrate at the alignment for the e fields itself. The superpotential for their alignment reads 

W e = S ei {4 - Ml) + 5 e . (j e) - M^) , (C.l) 

where i = 1, . . . , 11 and j = 12, 13. Note that for the sake of readibility we do not include any 
couplings. The driving fields S ti and S ej are total singlets so that terms like S € M^ are allowed. 
The .F-Term equations for the driving fields give, e.g. 

F Sei =el-Ml=0. (C.2) 

And since we assume that our fundamental theory is CP conserving the mass M ei is real (like 
the coupling parameters which are not shown) and hence the vev of e\ is real and non- vanishing. 
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For ei2 and ei3 this has to be slightly modified. For them we find three possible solutions, two 
of them complex and only one real. But we assume that the real solution is picked up, which 
could be preferred by higher order corrections, supergravity corrections or some low-energy soft 
terms in the scalar potential. To discuss this corrections in detail is beyond the scope of the 
current paper. 

Note also that all the S ei driving fields have the same quantum numbers and hence can mix 
with each other. In other words each of these driving fields could couple to each e field. We 
have chosen here the basis in which the superpotential has the above structure, which makes 
the alignment clear (see also the appendix of |42j). 

The same method can be applied to the real triplet and singlet flavons of our model, after 
we have fixed their alignment by some different kind of operators. But for the complex doublets 
(2', 2") and singlets (1', 1") we have to use other relations, because the representation squared 
cannot form a total singlet. 

Before we come to this complex representations we discuss the alignment for the flavons 
appearing in the neutrino sector (£, p, p) where this complication is absent 0. The superpotential 
for these flavons reads 

We,™ = x ^ eg + ipe9 + ^ + 5 « " M « ) + sp ( p2 + p 2 ~ M ') • (c - 3) 

The first thing to note here, is that we used the auxiliary flavon eg in the first set of operators 
involving the triplet driving field D^. Since eg appears in all three operators, it drops out in the 
F-term conditions, but nevertheless it is real and hence would just modify the value of the vev 
without introducing any phase. The F-term conditions are 



2£i 2 -2£ 2 £3 + £i(p + p) = 0, (C.4) 



3% 



dD 



£2 



2£2-2£i£3 + 60> + /5) = 0, (C.5) 



= 2f 3 - 266 + Up + P) = o , (C.6) 
W ^ -e i+ - Ml = , (C.7) 



^— - a -r *cks - 

p 2 + p 2 — M'p = . (Ci 



Besides the trivial solution £j = 0, i = 1,2,3, we find for the first three of these equations by 
cyclic permutations in £j the desired solution for which £j = £o 7^ if po = — Po- The fact that 
the vevs are non- vanishing and real can then be read off from the last two equations for which 
we used the method from [32] discussed above. 

Now we come to the most complicated part of the flavon alignment sector, the flavons present 
in the quark and charged lepton sectors. Although we have two different set of flavons, one for 
the up-quark sector at the one hand and one for the down-type quark and charged lepton sector 
at the other hand, we cannot separate their alignments completely. In fact, we found that the 
alignment is itself independent from each other but the simplest solution which we found to 
make all vevs real involves cross couplings between the two sectors. The flavon superpotential 



The alignment for the triplets follows the discussion in the seminal paper |43| . 
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reads 



w / = + ^" e + WC" + U - m~~C') + s ( (C'C" - 6 3 ) (C.9) 



+ ^ to - c'C") + ^ (#6 4 + 0C"e 5 ) + ^ (V') 2 66 + 0C'e 7 ) (C.10) 



+ 5^ (^V - M|) + ^ (V'C" + # - M ? <' + |c') (C.ll) 
+ 5 1 (^>"|-Ml 1 )+^((C') 2 6i 2 -(C / ) 2 6i3) , (C12) 

where in the last equations the cross couplings between the two sectors are written. We do not 
want to discuss here all the details of the alignment in detail, instead we will only discuss the 
phases of the vevs of the complex fields in a bit more detail. Nevertheless, we quote all the 
F-term conditions, in which it is then quite easy to plug in the flavon vevs from eqs. (|2.ip - (l2.3p 
and see that they form a viable solution. The F-term conditions read for the up sector 



for the down sector 



ei(20 2 -20 2 3 ) + e 2 2 C" = O, 
ei (2$- 20x0 3 ) + e 2 0iC" = O, 
ei(20!-2<M 2 ) + e 2 ^3C" = O, 
||f = (C") 2 - M ? ~,C' + % + 20!0 2 = 



dW f 
81)^ 
dWf 
d£>^ 



dWf 
dWf 



C'C" - M f -e 3 = , 

= 0, 

= e 4 (20 2 - 20 2 3 ) + e 5 02C" = 
= e 4 (202 -20 1 3 ) + e 5 0iC" = 
= e 4 (20|-20 1 02) + e 5 03C" = 
= 6 6 ((^ / ) 2 + 67^C' = 0, 



dW f 
dW f 

|^ = i6 6 («) 2 + e7 2 C / = O 



(1 - iJeoVW + er^iC = 
^ - ^» - Ml = , 
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9Wf = (c „,2 



(m c , + | V = , (C.26) 



and for the cross couplings between the two sectors 



dSi 

dWf 



- i>$f) € f-M 2 Sl =0, (C.27) 
(C') 2 ei2-(C') 2 ei3 = 0. (C.28) 



So how do we make the vevs of the complex representations real? Exemplary we discuss the 
complex singlets £" and From eqs. ()C.16p and ()C.17p we find a polynomial in £" 

(C") 3 + C"($ + 2&&) - M^M c -e 3 = , (C.29) 

which has a real solution (at least for a certain choice of parameters and plugging in the real vev 
of which we pick here. Then we know that (C") 3 ^ s rea l> while £' has the opposite phase of 
Q" so it is real as well. From eq. (]C.28p we then find £' to be real and from eq. ()C.26j) we obtain 
Q" to be real and all the singlet vevs are real. For the doublets a similar mechanism applies. 

The last alignment we want to discuss here is strictly speaking not an alignment. But since 
we have used adjoints of SU(5) in our operators to get the desired Yukawa coupling relations 
between the charged leptons and the down-type quarks we add here a mechanism which generates 
the vev of these adjoints and also show explicitly that they are real. For the fields and H'-j^ 
we can write down the following superpotential using the two driving fields S'^ and S'^ 

W 24 = S!U {H^H>> 4 - £ 2 ) + X (B&H&ew - £ 2 e n ) • (C.30) 

We see that the vev of £ triggers a vev for the two adjoint fields and even more these two vevs 
are directly related to the T' symmetry breaking scale. That means that in our model the GUT 
scale and the scale of T' coincide (up to some order one coefficients). In principle, we can again 
choose here between two different vevs for the adjoints. One pointing into the SM direction 
and the other one pointing into the SU{A) x U(l) direction and we assume the first option to 
be realised. We also note here, that the solution of the Doublet-Triplet-Splitting problem and 
hence the construction of the whole Higgs sector is clearly beyond the scope of this paper. 
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